Current rectification by simple molecular quantum dots: an ab-initio study 
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We calculate a current rectification by molecules containing a conjugated molecular group sand- 
wiched between two saturated (insulating) molecular groups of different length (molecular quantum 
dot) using an ab-initio non-equilibrium Green's function method. In particular, we study S(CH2)m- 
CioH6-(CH2)nS dithiol with Naphthalene as a conjugated central group. The rectification current 
ratio ~ 35 has been observed at m = 2 and n — 10, due to resonant tunneling through the molecular 
orbital (MO) closest to the electrode Fermi level (lowest unoccupied MO in the present case). The 
rectification is limited by interference of other conducting orbitals, but can be improved by e.g. 
adding an electron withdrawing group to the naphthalene. 

PACS numbers: 85.65.-Fh 



I. INTRODUCTION 

Effective current rectifier is a necessary element of elec- 
tronics circuitry. Aviram and Ratner suggested in 1974 
that in a molecule containing donor (D) and acceptor 
(A) groups separated by a saturated cr-bridge (insulator) 
group, the (inelastic) electron transfer will be more favor- 
able from A to D, rather than in the opposite direction 
[1] . It was noted more recently that the electron transfer 
from D to A should start at smaller bias voltage because 
the highest occupied molecular orbital (HOMO) is cen- 
tered on D group, whereas the lowest unoccupied molec- 
ular orbital (LUMO) is located on the acceptor group 
A. The onset of resonant tunneling in this case corre- 
sponds to the alignment of HOMO on the D group with 
LUMO on the A group under external bias voltage, which 
are closest in energy [2]. The molecular rectifiers (MR) 
synthesized in the laboratory were of somewhat differ- 
ent D — IT — A type, i.e. the "bridge" group was conju- 
gated [3,4]. The molecule, C16H33 — 7Q3CNQ, can be 
viewed as the fused naphthalene and TCNQ molecules 
with attached C16H32 alkane "tail" needed to provide 
enough van-der-Waals attraction between the molecules 
to assemble them into Langmuir-Blodgett film on a wa- 
ter. Although the molecule did show a rectification (with 
considerable hysteresis), it performed not like a conduc- 
tor but rather like an anisotropic insulator because of 
that large alkane "tail" . Indeed, the currents reported in 
[4] were on the order of 10~^^ A/molecule, which make its 
application impractical. It was recently realized that in 
this molecule the resonance does not come from the align- 
ment of HOMO and LUMO, since they cannot be decou- 
pled through the conjugated tt— bridge, but rather due 
to an asymmetric voltage drop across the molecule where 
HOMO and LUMO are asymmetrically positioned with 
respect to the Fermi level of the electrodes [5]. It is not 
clear if this mechanism has been observed. For instance. 



the contact between the alkane chain and electrode was 
certainly poor, via weak van-der-Waals forces. No tem- 
perature dependence of the resistivity was reported, but 
its large value may suggest that inelastic tunneling pro- 
cesses might have been involved. 

It is clear that the rectifying molecular films built with 
the use of LB technique relying on use of long aliphatic 
(or other insulating) "tails" will produce hugely resis- 
tive molecules, not suitable for moletronics applications. 
There are reports of rectifying behavior in other classes 
of molecules, e.g. molecules chemisorbed on electrodes, 
where an observed asymmetry of I-V curves is likely due 
to asymmetric contacts with electrodes [6,7], or asymme- 
try of the molecule itself [8] . 

To achieve a good rectification by a molecule with 
a reasonable conductivity, one should avoid using 
molecules with long saturated (insulating) groups, which 
would make them prohibitively resistive. Thus, a sug- 
gestion to use relatively short molecules with certain end 
groups to allow their self-assembly on a metallic elec- 
trode's surface seems to be very attractive [9]. In that 
previous paper the transport through a phenyl ring con- 
nected to electrodes by asymmetric alkane chains (CH2)„ 
has been calculated with the use of a semi-empirical 
tight-binding method. Since the microscopic parameters 
of the model are poorly known, especially hopping in- 
tegrals between molecule and electrodes, electrode work 
function and the affinity of the molecule (which may be 
strongly affected by bonding to a metal), the calcula- 
tions have been performed for a number of these param- 
eters. The calculations indicated a rectification ratio for 
-S-(CH2)2-C6H4-(CH2)„-S- molecule of about 100 for 
n = 10 with the resistance R « ISGil. It is difficult 
to know how quantitative these numbers are, because 
of a limited basis set and uncertainty in parameters of 
the model. For instance, an account for larger basis set 
for wave functions may lead to an increased estimate of 
the tunnel current, since the system would have more 



1 



channels to conduct current. Experimentally, it may 
be easier to synthesize somewhat different compounds 
with a naphthalene group as a "molecular quantum dot" , 
like -S-(CH2)„-CioH6-(CH2)„-S- [10], which we will 
also refer to as -S-(CH2)„-Naph-(CH2)„-S-. To ob- 
tain an accurate description of transport in this case, 
we employ an ab-initio non-equilibrium Green's function 
method [11]. 

The calculation shows that indeed the current rectifi- 
cation /-!_ //_ ~ 100 may be possible for some designs like 
-S-(CH2)3-Naph-(CH2)io-S-, where /+(_) is the for- 
ward (reverse) current. The difference between the for- 
ward and reverse voltages is, however, limited by other 
orbitals intervening into the conduction process. One 
needs the conducting orbital to be much closer in energy 
to the electrode Fermi level than the other ones (e.g. 
LUMO versus HOMO) and this energy asymmetry can 
be manipulated by "doping" the conjugated conducting 
part by attaching a donor (or acceptor) group, as will be 
shown below. 

The present calculation takes into account only elas- 
tic tunneling processes. Inelastic processes may substan- 
tially modify results in the case of strong interaction of 
electrons with molecular vibrations, see [15]. There are 
indications in the literature that the carrier might be 
trapped in a polaron state in saturated molecules slightly 
longer than we consider in the present paper [16]. 



II. RECTIFICATION BY MOLECULAR 
QUANTUM DOT 

The structure of the present molecular rectifier is 
shown in Fig. 1. The molecule consists of a central 
conjugated part (naphthalene) isolated from the elec- 
trodes by two insulating "arms" built from saturated 
aliphatic chains (CH2),i with lengths Li {L2) for the left 
(right) chain. The principle of molecular rectification by 
a molecular quantum dot is illustrated in Fig. 2, where 
the electrically "active" molecular orbital, localized on 
the middle conjugated part, is LUMO, which lies at the 
energy A above the electrode Fermi level at zero bias. 
The position of the LUMO, for instance, is determined 
by the work function of the metal g0 and the affinity 
of the molecule gX: ^ = Alumo — 9 (</" ^ x) ■ The po- 
sition of HOMO is given by Ahomo = Alumo - Eg, 
where Eg is the HOMO-LUMO gap. If this orbital is 
considerably closer to the electrode Fermi level Ep, then 
it will be brought into resonance with Ep prior to other 
orbitals when the molecule is biased in either forward or 
reverse direction (Fig. 2). It is easy to estimate the for- 
ward and reverse bias voltage assuming that the voltage 
mainly falls off at the saturated (wide band gap) parts of 
the molecule with the lengths Li (L2) for the left (right) 
barrier. Fig. 2. Since the voltage drops proportionally to 
the lengths of the barriers, 14(2) cx: -^1(2)7 we obtain for 




FIG. 1. Schematic diagram of typical Au-molecular recti- 
fier- Au device structure analyzed within our ab-initio formal- 
ism. Electrodes are constructed as the 111 orientation, with 
the terminal sulfurs of the molecule occupying the triangular 
hollow site. Li (L2) indicates the insulating barrier length 
separating the molecular conducting unit from the left (right) 
electrode. 



the partial voltage drops Vn2) = ^^1(2) /^i where V is 
the bias voltage. Forward and reverse voltages are found 
from the resonance condition, which gives 



Vf = -(1 + 0. 

Vp/Vr = i = L1/L2 



(1) 

(2) 
(3) 



where q is the elementary charge. A significant difference 
between forward and reverse currents should be observed 
in the voltage range Vp < \V\ < Vr. For our systems 
in question the parameter A = 1.2 — 1.5 eV, Fig. 3, 
the value determined by the work function of the elec- 
trodes and electron affinity of the conjugated part of the 
molecule, as well as by the strength of a surface dipole 
layer. For unlike electrodes with the work functions q(f>i 
and q(f>2 there will be a built-in potential in the system 
qcfib — q {4>2 — 0i), Fig. 2(d). In this case the expressions 
for the forward and reverse bias voltage are modified as 
follows: 



Vp 



+ 0b (1 + 0, 



A / 1 

Vfl = 0b + — 1 + - 
9 V ^ 



(4) 
(5) 



Thus, the operating bias "window" in this case 
Vr-Vp = {Vr - Vf)o - 



(6) 



Therefore, it slightly (^ < 1) increases when (jjh < {(jji > 
02, left work function is larger), and shrinks down when 
0fc > (right work function is larger). One would like to 
use as asymmetric a molecule as possible to increase the 
operating voltage, and in this case the difference of the 
work functions of the electrodes progressively becomes 
less important. 
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FIG. 2. Schematic band diagram of the molecular rectifier. 
Middle narrow-band part represents the naphthalene conju- 
gated group. The barrier height Uo — 4.9 eV (present theory), 
4.8eV (experiment [16]). (a)-(c) Unbiased, forward, and re- 
verse biased rectifier with similar electrodes, (d) Unbiased 
rectifier with dissimilar electrodes. At equilibrium, there is a 
voltage drop gVi(2) across the left (right) barrier. 



The main transport features of the present system can 
be qualitatively explained within a model of resonant 
tunneling through a localized state (molecular orbital) 
in the barrier. The current can be written as 



2q 
h 



I dET{E) [f{E) - f{E - 



qV)], 



(7) 



where the transmission T[E) at zero temperature is given 
by (for two spin directions) 



T{E) 



Li R 



(8) 



where Emo ^ Ef + A - qVi = Ep + A - qVLi/L is 
the energy of the molecular orbital, relative to the Fermi 
level of the left electrode. The transmission is maximal 
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FIG. 3. Transmission coefficient (top) and local den- 
sity of states (bottom) versus energy E for rectifiers 
-S-(CH2)2-CioH6-(CH2)„-S-, n = 2,4,6,10. A indicates 
the distance of the closest MO to the electrode Fermi energy 
{Ef = 0). 





FIG. 4. Projected density of states N{E, r) , onto a 2-D grid 
defined by the plane of our central molecule, for the naphtha- 
lene series, n = 2, 4, 6, 10 from left to right, respectively. As 
the alkane chain increases in length, the coupling of the naph- 
thalene group to the right electrode decreases. 



and equals 1 when E = Emo and = F^j, which corre- 
sponds to a symmetric position of the central conjugated 
part with respect to the electrodes. We can estimate 
^L{R) ~ /D = Foe~2''^i(2)^ where t is the overlap inte- 
gral between the MO and the electrode, D is the electron 
band width in the electrodes, k the inverse decay length 
of the resonant MO into the barrier. The latter quantity 
can be estimated from Tt^k^ /2m* — Uq, where Uq is the 
barrier height 4.8 eV in alkane chains (CH2)„ [16], 
« 5 eV from DFT calculations], and m* < 1 the effective 
tunneling mass. The current above the threshold (when 
the resonant tunneling condition is met) is given by 



2q 



Toe 



-2kL2 



(9) 



We see immediately that increasing the spatial asymme- 
try of the design (i.e. increasing L2/L1) changes the op- 
erating voltage range hnearly, Vr — Vp « {A/q) L2/L1, 
and brings about an exponential decrease in current [9]. 
This severely limits the ability to increase the rectifica- 
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tion ratio while simultaneously keeping the resistance of 
the molecule at a reasonable value. One therefore needs 
a quantitative analysis of all microscopic details, includ- 
ing molecule-electrode contact, and potential and charge 
distribution, in order to make comparisons with the data. 

III. METHOD 

Since we study relatively short molecules, electron 
transport is likely to proceed elastically. We shall treat 
the elastic tunneling processes from first principles in 
order to evaluate all the relevant microscopic charac- 
teristics within the density functional theory. This al- 
lows us to check the accuracy of semi-empirical meth- 
ods used previously, and find results for new objects 
of experimental interest, like the present molecular 
quantum dot with a naphthalene conjugated middle 
group. We use an ab-initio approach that combines the 
Keldysh non-equilibrium Green's function (NEGF) [12] 
with pseudopotential-based real space density functional 
theory (DFT) [11]. In contrast to standard DFT work, 
our NEGF-DFT approach considers a quantum mechani- 
cal system with an open boundary condition provided by 
semi-infinite electrodes under external bias voltage. The 
main advantages of our approach are (i) a proper treat- 
ment of the open boundary condition; (ii) a fully atom- 
istic treatment of the electrodes and (iii) a self-consistent 
calculation of the non-equilibrium charge density using 
NEGF. 

For the open devices studied here, our system is di- 
vided into three regions, a left and right electrode and a 
central scattering region, which contains a portion of our 
semi-infinite electrodes (see Fig. 1). By including atomic 
layers of the electrode in our scattering region, we provide 
a proper treatment of the electrode-molecule interaction, 
and allow for the effective screening of the molecule away 
from our molecular junction. External bias voltage can 
then bo applied as a boundary condition for the Poisson 
equation in our central scattering region. The Hartree 
potential in the scattering region is solved using a multi- 
grid real-space numerical approach. Finally, the non- 
equilibrium density matrix p is constructed via NEGF 
[12], 

p=^.JdETrG<{E), (10) 

where the lesser Green's function is determined from the 
operator relation 

G<=G^Y,<G^, (11) 

via the retarded (advanced) Green's functions G^^^^ and 
T,^ the self-energy part describing an injection of charge 
into the scattering region (molecule) from the electrodes. 
This process is characterized by the usual open channel 



representation of the scattering states with momenta fc" 
(A;") in the left (right) electrodes, where index n enu- 
merates all open channels (available Bloch states of the 
electrons in the electrodes). The Fermi functions f{k") 
and /(fc") define the occupied states in the leads to be 
accounted for in evaluation of the electron transport. To 
evaluate the Green's function, one treats the disconnected 
system consisting of the left (right) lead marked by index 
I (r) and the scatterer c. The transport Green's function 
is then found from the Dyson equation 

{G^)-' = {G^)-'-V, (12) 

where the unperturbed retarded Green's function is de- 
fined in operator form as 

(G^)"' = {E + iO)S-H, (13) 

/ Hi,i Hi,, \ 
H = H,,i H,,r ■ (14) 

\ Hr^c Hr^r / 

H is the Hamiltonian matrix, Hc,i = ^, Hr.c = H}. S 
is the overlap matrix, Sij = {xi\Xj) ^oi non-orthogonal 
basis set orbitals Xi; and the coupling of the scatterer 
to the leads is given by the Hamiltonian matrix V = 
diag[S;^(, 0, Sr-,r]- In all formulas above it is assumed 
that matrix elements are given in terms of the non- 
orthogonal basis set. 

The self-energy part S<, which is used to construct the 
non-equilibrium electron density in the scattering region, 
is found from 

S< = -2i Im [/(^)Sy + f{E + qV)i:r,r] ■ (15) 

This expression accounts for the steady charge "flowing 
in" from the electrodes. The transmission probability, 
which determines the current according to Eq. (7), is 
given by 

T{E) = 4Tr [Im (S,,,) Gf, Im (S,,^) G^^i] . (16) 

The local density of states, which is a very important 
quantity to characterize the spatial character of the wave 
functions, can be obtained from 

N{E,f) = ^Tr [G<{E) \x{r)) (x(r')|] . (17) 

To construct the system Hamiltonian, we define the 
atomic cores using standard norm- conserving nonlo- 
cal pseudopotentials [13], and expand the Kohn-Sham 
wave functions in an s-,p-,d- real space atomic orbital 
basis [11,14]. External bias is incorporated into the 
Hartree potential, and in this way the non-equilibrium 
charge density is iterated to self-consistency, after which 
the current- voltage (I-V) characteristics have been com- 
puted. 
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FIG. 5. Transmission through molecular 

rectifiers -S-(CH2)™-CioH6-(CH2)6-S-, m = 0, 1, 2. The 
peak corresponds to transmission through the LUMO state 
localized on the naphthalene ring. 

IV. RESULTS 

The physical picture of the transmission through a 
series of rectifiers — S-(CH2)m-CioH6 — (CH2)„— S— for 
m = 2 and n — 2,4,6,10 is illustrated in Figs. 3,4. 
The transmission coefficient T{E) (top panel) and the 
density of states N{E) = / dVN{E, f) (bottom panel) 
are shown at zero bias voltage in Fig. 3. The Fermi 
level in the electrodes is taken as the energy origin. It 
follows from this figure that indeed the LUMO is the 
closest molecular orbital transparent to electron trans- 
port, and it is higher in energy than Ep by an amount 
A = 1.2 — 1.5 eV depending on the molecule. Since the 
molecule is asymmetric, even at the resonance energy we 
have T{Elumo) 0-5j Fig. 3. This can be traced back 
to the Breit-Wigner formula (8), which for Ffl < Tl gives 
T{Elumo) ~ 4Ffl/FL < 1, and it falls off as the asym- 
metry increases. There the transmission through nominal 
HOMO state is neghgible, but HOMO-2 state conducts 
very well. As we shall see below, the HOMO-2 defines the 
threshold reverse voltage Vr, thus limiting the operating 
voltage range (6). 

The maps of the densities of states N{EmMO,^) in 
Fig. 4 show the development of the tunneling gap be- 
tween the naphthalene group and more distant electrode 
separated by alkane chain (CH2)„, n = 2 — 10. The 
symmetry and weight distribution of LUMO in the de- 
vice is very similar to LUMO in the isolated molecule, 
telling us that the character of the molecular wave func- 
tions is preserved to a high degree. Obviously, in con- 
tact with electrodes molecular orbitals acquire some finite 
width r = Fl -t-F/j, and for highly asymmetric molecules 
F « Fi, so it is mainly defined by the transparency of 
the shorter (left) insulating group (CH2)m. One needs a 
rather small width F which defines the sharpness of the 
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FIG. 6. Potential drop as a function of position along the 
length (transport direction) of our molecular device, for the 
n = 2, 6, 10 molecular rectifiers, at an applied bias voltage of 
IV. The positions of the naphthalene ring inside the junction 
are indicated by the stick figures, and the respective insulat- 
ing barrier lengths are shown as horizontal lines. Curves are 
vertically shifted for clarity. 

voltage threshold for the current turn-on [9]. We have 
calculated the molecules with m = 0, 1, 2 with the results 
for F,„ = 130, 60, and 12 meV (Fig. 5). The last value for 
the left insulating group (CH2)2 seems to be most rea- 
sonable, since the resonant peak is rather narrow yet the 
molecule remains conductive. However, in some cases the 
short insulating group (CH2)3 may be preferential. This 
conclusion is similar to the one reached on the basis of the 
tight-binding model, but the peak positions and size are 
substantially different in both cases, the peak size being 
much larger in the case of the present DFT calculations, 
Fig. 5. 

Our assumption that the voltage drop is proportional 
to the lengths of the alkane groups on both sides, is quan- 
tified by the calculated potential ramp in Fig. 6. It is 
close to a linear slope along the (CH2)„ chains. There is 
also a smaller but still noticeable voltage drop across the 
naphthalene group, which can be viewed as a fragment 
of wide-band semiconductor. 

The threshold bias voltages can be obtained from 
Fig. 7, showing the evolution of the LUMO versus 
bias voltage with respect to the chemical potentials 
of the right (left) electrodes fiR(L) = Ep + {—)qV/2. 
The forward voltage corresponds to the crossing of the 
LUMO(y) and fiR{V), which happens at about 2V. Ex- 
trapolated crossing of LUMO and /zl is at large nega- 
tive voltages but, unfortunately, this large desirable dif- 
ference between Vp and Vr does not materialize. Al- 
though the LUMO defines the forward threshold voltages 
in all the molecules studied here, the reverse voltage is 
defined by the HOMO-2 for "right" barriers (CH2)„ with 
n — 6, 10. The I-V curves are plotted in Fig. 8 and the 
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FIG. 7. Evolution of the rectifier LUMO energy state as a 
function of applied bias voltage. Solid black lines represent 
the electrode chemical potentials, which are symmetrically 
shifted by the applied voltage j-iR(L) ~ Ef + {—)qV/2. The 
forward voltage threshold corresponds to the crossing of the 
LUMO energy with the right electrode chemical potential. 
In the case of CN-doped naphthalene ring both forward and 
reverse threshold voltages correspond to LUMO crossing the 
electrode Fermi level, as in the ideal case. Fig. 2. In other 
systems the HOMO defines the reverse voltage threshold, thus 
reducing the operating voltage range. 

TABLE I. The parameters of a set of the naphtha- 
lene-based rectifiers with various right insulating groups 
(CH2)n, Fig. 1. I+/I- the currect rectification ratio, R the 
molecule resistance at bias voltage V — 2.5V. 



n 


A (eV) 


Vf (V) 


Vr (V) 


I+/I- 


R (un) 


2 


1.23 


2.45 


2.45 


1 


2.20 


4 


1.37 


2.50 


3.00 


5 


20.1 


6 


1.39 


2.15 


3.13 


15 


90.9 


10 


1.50 


2.10 


2.85 


35 


2780 



corresponding parameters are listed in the Table. We see 
that the rectification ratio for current in the operation 
window /+ / 1- reaches the maximum value of 35 for "2- 
10" molecule {m = 2, n — 10). It was estimated to be 
about 100 for "2-10" molecule with a phenyl ring as a 
conjugated central part [9]. We have checked a series of 
molecules with central phenyl ring and, unfortunately, do 
not find any significant rectification in this case, because 
the Fermi level in DFT calculations appears to lie close 
to mid-gap. Thus, the tight-binding results for large rec- 
tification in one-phenyl ring molecules might have been 
the results of the model. The estimated resistance of the 
molecule in DFT is substantially less (2.8 GO) compared 
to the tight-binding model with a phenyl ring as a con- 
jugated central unit (13 GO), see Table. 

One can manipulate the system in order to increase 
the energy asymmetry of conducting orbitals (reduce A). 
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FIG. 8. I-V curves for naphthalene recti- 
fiers -S-(CH2)2-CioH6-(CH2)„-S-, n = 2,4,6,10. The 
short-dash-dot curve corresponds to a cyano-doped (added 
group -C=N) n = 6 rectifier. 

To shift the LUMO towards Ep, one needs in effect to 
gate the naphthalene group by positive voltage. A chem- 
ical way of doing this would be to attach some chemical 
group, like -C=N in the present case. This group with- 
draws a fraction of an electron and acquires the negative 
charge —d = — O.llg. The charge is donated mainly by 
the naphthalene group, which acquires a positive charge 
and orbitals on the molecule shift to lower energies al- 
most rigidly by about 0.6 eV, Fig. 9. On the resulting I-V 
curve, both thresholds Vf(_r) now correspond to LUMO 
crossing cither the right or left electrode Fermi levels (see 
Fig. 7), and the operating voltage range is modified sig- 
nificantly. 

V. CONCLUSIONS 

We have presented parameter-free DFT calculations 
of a class of molecular quantum dots showing a cur- 
rent rectification at a reasonable current density through 
the molecule. The explored mechanism of rectification 
is different from Aviram-Ratner diodes [1] and relies on 
considerable asymmetry of the spatial composition and 
energy structure of conducting orbitals of the molecule 
[9]. Since there is no empirical parameters, the predic- 
tions are quantitative and can be tested in a control- 
lable fashion [10]. We have studied the molecular quan- 
tum dots with the naphthalene central conjugated group, 
and found a rectification of about 35 for — S-(CH2)2- 
CloH6— (CH2)6— S— molecules. The modest rectification 
ratio is a consequence of the resonant processes at these 
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FIG. 9. Left: Energy level positions for undoped and 
cyano doped (-C=N) symmetric (m,n = 2) naphthalene rec- 
tifier. Right: Transmission T{E) for the doped (solid) and 
undoped (dotted) rectifiers at zero applied bias voltage. The 
-C=N substituent shifts the electronic levels near Ep to lower 
energy, increasing the energy asymmetry and hence decreas- 
ing A. 



molecular sizes. The present study suggests also that the 
rectification ratio in the case of a central single phenyl 
ring, studied previously [9], is dramatically overestimated 
by the tight-binding model. 

The rectification ratio is not great by any means, but 
one should bear in mind that this is a device necessarily 
operating in a ballistic quantum-mechanical regime, be- 
cause of small size. This is very different from present 
Si devices with carriers diffusing through the system. 
As they become smaller, however, the same effects, as 
those discussed here for molecular rectifiers, will eventu- 
ally take over and will tend to diminish the rectification 
ratio. In this regard the present results are indicative of 
the problems which will ultimately be faced by Si devices 
of molecular proportions. 

It is important that the parameters of the device, 
like the threshold voltages for forward and reverse bias, 
can be significantly modified by "chemical" doping. In 
the present case adding the -C=N electron withdraw- 
ing group has scaled the threshold range by a factor of 
about two. There may be other modifications that could 
improve/add functionality of the molecular devices by 
varying central, end, and side groups of the molecule. 
The present results illustrate very clearly that the opti- 
mization of the parameters of the molecular devices is 
typically facing the problem that some parameters, like 



level widths, vary exponentially, while others, like the 
voltage thresholds, vary as a power law- as a function of 
the molecular design. Although the approximate band 
diagrams illustrating the level structure. Fig. 2, are very 
helpful for qualitative description, the atomistic details 
of the electron density distribution are extremely impor- 
tant, Fig. 4. This necessitates the use of ab-initio meth- 
ods, without which one may obtain qualitatively wrong 
results. 

We have found that in the present molecules signifi- 
cant rectification may be obtained in relatively short self- 
assembled molecules. Self-assembly provides for good 
contact with electrode and relatively small resistance 
per molecule. In existing studies of Aviram-Ratner-type 
molecules they were assembled by Langmuir-Blodgett 
deposition and are prohibitively resistive by design be- 
cause of a long aliphatic tail needed for the assembly 
[3,4]. Thus, for apphcations the self-assembly of shorter 
molecules on metallic electrodes seems to be essential. 
One should keep in mind that the molecular rectification, 
here estimated for ideal contacts at zero temperature, will 
be reduced at finite temperature and due to inevitable 
disorder in the contact area and in the molecular film 
[17]. Indeed, the optimal width of the level was found to 
be about 12 meV, which is just a half of the room tem- 
perature, so the temperature (and static disorder) will 
significantly broaden the level and reduce the rectifica- 
tion. In terms of simplest cross-point memories one may 
use a large number of molecules on a pitch. This may al- 
low for use of longer, better rectifying molecules. There, 
however, the inelastic processes will play a progressively 
more significant role [16] and will further reduce the recti- 
fication. Further studies of molecular rectifiers should be 
focused on these outstanding issues in conjunction with 
experimental studies. 
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